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Abstract 

This paper is the widely extended version of the pubhcation, appeared in Proceedings of IS- 
SAC'2009 conference (Andres, Levandovskyy, and Martin- Morales, 2009). We discuss more de- 
tails on proofs, present new algorithms and examples. We present a general algorithm for com- 
puting an intersection of a left ideal of an associative algebra over a field with a subalgebra, 
generated by a single element. We show applications of this algorithm in different algebraic sit- 
uations and describe our implementation in Singular. Among other, we use this algorithm in 
computational D-module theory for computing e. g. the Bernstein-Sato polynomial of a single 
polynomial with several approaches. We also present a new method, having no analogues yet, 
for the computation of the Bernstein-Sato polynomial of an affine variety. Also, we provide a 
new proof of the algorithm by Briangon-Maisonobe for the computation of the s-parametric an- 
nihilator of a polynomial. Moreover, we present new methods for the latter computation as well 
as optimized algorithms for the computation of Bernstein-Sato polynomial in various settings. 



1. Introduction 

This paper extends Andres, Levandovskyy, and Martin-Morales (2009) by many de- 
tails, proofs, algorithms and examples. In this paper we continue reporting (Levan- 
dovskyy, 2006; Levandovskyy and Martin-Morales, 2008; Andres et al., 2009) on our 
advances in constructive Z?-module theory both in theoretical direction and also in the 
implementation, which we create in Singular. 

Our work on the implementation of procedures for _D-modules started in 2003, mo- 
tivated among other factors by challenging elimination problems in non-commutative 
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algebras, which appear e. g. in algorithms for computation of Bernstein-Sato polyno- 
mials. We reported on solving several challenges in Levandovskyy and Martin-Morales 
(2008). A non-commutative subsystem Singular:Plural (Greuel et al., 2006) of the 
computer algebra system SINGULAR provides a user with possibilities to compute nu- 
merous Grobner bases-based procedures in a wide class of non-commutative G-algebras 
(Levandovskyy and Schonemann, 2003). It is natural to use this functionality in the 
context of computational D-modulc theory. 

As of today, the D-module suite in Singular consists of three libraries: dmod.lib, 
dmodapp.lib and bfun.lib. Moreover, gmssing.lib (Schulze, 2004) contains some so- 
phisticated (and hence fast) and useful procedures, e. g. bernstein for the computation 
of the local Bernstein- Sato polynomial of an isolated singularity at the origin. There 
are many useful and flexible procedures for various aspects of D-module theory. These 
libraries are freely distributed together with Singular (Greuel et al., 2005) since the ver- 
sion 3-1-0, which was released in April 2009. More libraries are currently under develop- 
ment, among them procedures for computing the restriction, integration and localization 
of D-modules. 

There are several implementations of algorithms for D-modules, namely the exper- 
imental program KAN/smI by N. Takayama (Takayama, 2003), the bfct package in 
Risa/Asir (Noro et al, 2006) by M. Noro (Noro, 2002) and the package Dmodules.m2 
in MACAULAY2 by A. Leykin and H. Tsai (Tsai and Leykin, 2006). To the best of our 
knowledge, there is ongoing work by the CoCoATcam (2009) to develop some D-module 
functionality as well. We aim at creating a _D-module suite, which will combine flexibil- 
ity and rich functionality with high performance, being able to treat more complicated 
examples. 

We continue comparing our implementation (cf. Section 7.1) with the ones in the 
systems AsiR and Macaulay2, see Levandovskyy and Martm-Morales (2008) for earlier 
results. 

In this paper, we address the following computational problems: 

• s-parametric annihilator of / € lK[a;i, . . . , x„] 

• Bernstein- Sato ideals for f = fi ■ . . . ■ /,„ 

• 6-function with respect to weights for an ideal in D 

• Global Bernstein-Sato polynomial of / 

• Bernstein- Sato polynomial for a variety 

In Section 3.3, we give a new proof for the algorithm by Briangon-Maisonobe for 
computing Ann£)[s] announced in Levandovskyy and Martm-Morales (2008). More- 
over, using the same technique we design a new algorithm for the computation of the 
Bernstein-Sato polynomial of an affine variety, following the paper Budur et al. (2006), 
and prove its correctness. 

We develop the method of principal intersection 4.11 in the general context of IK- 
algebras and discuss its improvements. This algorithm is especially useful for problems 
of il'-module theory, since it allows to replace a generally hard elimination with Grobner 
bases by the search for a K-linear dependence of a sequence of normal forms. The al- 
gorithm is applied in Section 5 to two main methods for computing Bernstein-Sato 
polynomials as well as to solving 0-dimensional systems in commutative rings and to 
the computation of central characters in Section 4.2.1. Moreover, we describe a folklore 
method for computing Bernstein-Sato polynomial via annihilator (using, however, prin- 
cipal intersection instead of Grobner-based elimination) and prove in Lemma 6.4, that is 
is more efficient than the usual one. 
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The generalization of principal intersection approach to the case of more general sub- 
algebras we discuss in Section 4.3. 

1.1. Notations 

Throughout the article K stands for a field of characteristic zero. By R we denote the 
polynomial ring K.[xi, . . . , Xn] and by / € i? a non-constant polynomial. 

We consider the n-th Weyl algebra as the algebra of linear partial differential operators 
with polynomials coefficients. That is D„ = = K(xi, . . . , x„, 9i, . . . , 9„ [ {diXi = 

Xidi + 1, diXj = Xjdi, i ^ j}). We denote by -D„[s] = Dn{R) <8)k 1K[si, . . . , s„] and drop 
the index n depending on the context. 

The ring ii is a natural £'„(i?)-module with the action 

Xi • f{x) = Xi ■ f{x), di • f{x) = ^1^- 

Working with monomial orderings in elimination, we use the notation x ^ y for "x is 
greater than any power of y" . 

Given an associative K-algebra A and some monomial well-ordering on A, we denote 
by lm(/) (resp. lc(/)) the leading monomial (resp. the leading coefficient) oi f G A. 
Given a left Grobner basis G C A and f € A, we denote by NF(/, G) the normal form of 
/ with respect to the left ideal a{G). We also use the shorthand notation h -^h f (and 
/i — > /, if is clear from the context) for the reduction oih G Ato f G A with respect to 
the set H. If not specified, under ideal we mean left ideal in a IK-algebra, and by Grobner 
basis a left Grobner basis. For a, b in some K-algebra A, we use the Lie bracket notation 
[a,b] := ab — ba as well as skew Lie bracket notation [a,b]k := ab — k ■ ba fov k € K* . 

We say, that a proper subalgebra S of an associative K-algebra A is a principal sub- 
algebra, if there exists 5 € ^ \ K, such that S = K[g] . 

Let M be an ^-module, then we denote the Gel'fand-Kirillov dimension of M (see 
McConnell and Robson (2001) for the details and Bueso et al. (2003) for algorithms) 
by GK. dim(Af). Recall, that a module M is called holonomic, if GK. dim(A/ AnnM) = 
2 • GK. dim(Af). We prefer this general definition, since it concides with the classical way 
of defining holonomy in Weyl algebras and it is incomparably more general to the latter. 
In particular, armed with this general definition we can speak on holonomic modules 
over any G-algebra. 



2. Preliminaries 

It is convenient to treat the algebras we deal with in the bigger framework of G-algebras 
of Lie type. 

Definition 2.1. Let A be the quotient of the free associative algebra K(3;i, . . . ,.t„) by 
the two-sided ideal /, generated by the finite set {xjXi — XiXj — dij} VI < i < j < n, 
where dij G K[a;i, . . . , a;„]. The algebra A is called a G-algebra of Lie type (Levandovskyy 
and Schonemann, 2003), if 

• \/l<i<j<k<n the expression dijXk — Xkdij + xjdik — dikXj + djkXi — Xidjk 
reduces to zero modulo / and 

• there exists a monomial ordering -< on K[a;i, . . . , a;„], such that Im(dij) -< XiXj for 
each i < j. 
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These algebras were also studied in Kandri-Rody and Wcispfcnning (1990) and Bueso 
et al. (1998) by the names PBW algebras and algebras of solvable type. 

Recall the algorithm for computing the preimage of a left ideal under a homomorphism 
of G-algebras from Levandovskyy (2006). 

Theorem 2.2 (Preimage of a Left Ideal, Levandovskyy (2006)). 

Let A, B be G-algebras of Lie type, generated by {xi | 1 < « < n} and {yj | 1 < j < m} 
respectively, subject to finite sets of relations Ra, Rb as in Def. 2.1. Let (}> : A ^ B he a. 
homomorphism of K-algebras. Define to be the (A, ^)-bimodule A{{xi — 4>{xi) \ 1 < 
i < i^})a C A®k-B. Suppose, that there exists an elimination ordering for B on A^-^B, 
satisfying the following conditions 

1 < i < n, 1 < j < TO, \\n{yj(f){xi) - 4>{xi)yj) -< Xiyj. 

Then there are the following statements. 

1) Define A B to be the K-algebra generated by {xi, ... , a;„, j/i, . . . , ym} subject to 
the finite set of relations composed of Ra,Rb and {yjXi — Xiyj — yj(j){xi) + (j){xi)yj}. 
Then A (g)^ i? is a G-algebra of Lie type. 

2) Let J C B be a left ideal, then 

0-i(j) = (/^ + j)nA c A0^BnA 

Moreover, this computation can be done by means of elimination. 

The following proposition is a reformulation of Theorem 2.2, adopted to the situation, 
which is often encountered in context of £)-modules. 

Proposition 2.3. Let Ai,Bi,C be G-algebras of Lie type and if : Ai ^ Bi he a 
homomorphism of K-algebras. Consider the following data: 

A = C'SikAi, B = C'SikBi, (j)=lc'Siip:A^B, 
E = A<S)iB, E' = C^k{Ai^IBi). 

Then A c E' c E and for a left ideal J C B we have: 

(1) {EI^ + EJ) nE' = E'l^ + E'J. 

(2) (J) = {E'I^ + E'J)nA. 

Moreover, the second intersection can be computed using Grobner bases, provided there 
exists an elimination ordering for Bi on £" compatible with the G-algebra structure of 
E'. 

In the proofs we quite often use the following. 

Lemma 2.4 (Generalized Product Criterion, Levandovskyy and Schonemann (2003)). 
Let A be a G-algebra of Lie type and f,gGA. Suppose that lm(/) and lin{g) have no 
common factors, then spoly(/, 5) reduces to [f,g] with respect to the set {f,g}. 

3. s-paremietric Annihilator of / 

Recall Malgrange's construction ior f = fi ■ . . . ■ fp £ R — K[xi, . . . ,x„]. Let Dn = 
Dn{R), T = K[ti, ...,tp] and D'^ := Dp{T). Consider the (p+n)-th Weyl algebra £)p+„ = 
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-D„(8)K-Dp. Moreover, consider the following left ideal in Dp+n, called the Malgrange ideal 

P Qf. 

If-={{ tj - fj,di + -^^tj' 1 < J < 1 < « < ^ } )• 
j=i ''^^ 

Then for s = (si, . . . , Sp) we denote := V . . fp". Let us compute 

If nK[{tjdtj}]{xi,dxi I [dxi,Xi\ = 1) C Dn[{tjdtj}] C fp+n 

and furthermore, replace tjdtj with —sj — 1. The result is known (e. g. Saito et al. (2000)) 

to be exactly Ann^jg] C D[s]. 

There exist several methods for the computation of the s-parametric annihilator of 

3.1. Oaku and Takayama 

The algorithm of Oaku and Takayama (Oaku, 1997a, b,c; Saito et al., 2000) was de- 
veloped in a wider context and uses homogenization. With notations as above, let H := 
Dn ®K Dp K[ui, . . . , Up, vi,. . . ,Vp]. Moreover, let / below be the {u, u)-homogenized 
Malgrange ideal, that is the left ideal in H 

I = (^{tj - Ujfj,d, + ^ y^Ukdtj,UjVj - 1}^ . 

Oaku and Takayama proved, that Ann^i^js] (/*) can be obtained in two steps. At first 
{uj,Vj} are eliminated from / with the help of Grobner bases, thus yielding I' = I C\ 
{Dn <8iK D'p). Then, one calculates I' D (D„ ®^ K[{—tjdtj — 1}]) and substitutes every 
appearance of tjdtj by —Sj — 1 in the latter. 

3.2. Briangon and Maisonobe 

Consider Sp = K{{dtj,Sj} \ dtjSk — Skdtj — djkdtj) (the p-th shift algebra) and 
B = Dn (S)k Sp. Moreover, consider the following left ideal in B: 

i = ({sj + fjdtj,di + p^^mk}y 

Briangon and Maisonobe proved in Briancon and Maisonobe (2002) that Ann£)^[g](/*') = 
/ n Dn[s] and hence the latter can be computed via the left Grobner basis with respect 
to an elimination ordering for {dtj}. 

3.3. A new proof for Briangon- Maisonobe 

By using the Preimage Theorem 2.2 we give a new, completely computer-algebraic 
proof for the method of Briangon-Maisonobe 3.2. 

Let A := Dn[s] = K{{sj,Xi,di} | diXi = Xidi + 1), and 

B := K{{tj, dtj, Xi, di} \ {diXi = Xidi + 1, dtjtj = tjdtj + 1}). 

Thus in the notations of Proposition 2.3, C = Dn, A\ = K[s], Bi = K{{tj,dtj \ dtjtj = 
tjdtj + 1}) and A — C ® Ai,B — C ® Bi. Consider the algebraic Mellin transform (cf. 
Saito et al. (2000)) : Ai ^ Bi, Sj ^ -tjdtj - 1. 
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Hence = {{tjdtj + Sj + 1}) C Ai <E>k-^i —'■ ^' ■ Since \tk, Sj] = 6jktj and [dtk, sj] = 
—6jkdtj, the ordering conditions of Theorem 2.2 take the form tj -< Sjtj,dtj -< sjdtj, 
which are satisfied if and only if 1 < tj , dtj , Sj . 

By Proposition 2.3, for any L C B, (j)~^{L) = (/^ + L) Ci A. Hence, 

P Qf 

h + L= {{tj -fj,di + Y, ^dtj,tjdtj + Sj + 1}) 

P 

= {{tj - fj, di + J2 fj^tj + Sj}) 

because tjdtj + Sj + 1 reduces to 

tjdtj + Sj + 1 — dtj ■ {tj — fj) = fjdtj + Sj e + L. 

Lemma 3.1. Consider an ordering -<t, which satisfies the property {tj} ^ {xi}, 
{di,Sj} » {xi,dtj}. Moreover, set 

9i ■■=^^ + J2 jr^^tk, Si := {tj - fj,gi}, S2 := 5i U {sj + fjdtj} C E' . 
fe=i 

Then Si and S2 are left Grobner bases with respect to -<t- 

Proof. We run Buchberger's algorithm by hand. There are only three kinds of critical 
pairs we have to consider. Due to the ordering property, for each pair the generalized 
Product Criterion is applicable. Hence, we need to compute just the Lie brackets of 
members of pairs. 

1. [U- f^,tk- fk]=()■ 
2. For pairs {g^gk) computing \gi,gk\ yields 



Since [di, ^] = = [dk, |^], s^o\y{gi,gk) reduces to zero. 

3. For mixed pairs {tk — fk,gi) we have 

[tk - fk,g^] = ^ If^t^'^'^^j] - [/'^'^»] = 0- 

j * 

Hence, is a left Grobner basis. Now, in ^2 there are three new kinds of critical 

pairs to consider and for all of them we can apply the generalized Product Criterion. 

4. [tk - fk, Sj + fjdtj] = [tk: Sj] + fAtk,dtj] - [/fe, Sj] - [fk, fjdtj] = 6jk{tk - fj) ^ 0. 

5. [si + fidti, Sj + fjdtj] = fj[si, dtj] - fi[sj,dti] = 0. 
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6. Finally, 





= S9tj-[di,fj]dtj = 0. 



So, 52 is a left Grobner basis. □ 

We want to eliminate both {tj} and {dtj} from + L. As we see above, by using an 
elimination ordering for {tj} we proved above that ^2 is a Grobner basis. So, the elimi- 
nation ideal is generated by iSa := 52 \ {ti — fi}. Hence we can proceed with eliminating 
{dtj} from ^3, which is exactly the statement of the Briangon-Maisonobe algorithm in 
Section 3.2. 

3.4- Bernstein-Sato ideals for f = fi ■ ■ ■ ■ ■ fm 

Comparing the effectiveness of the algorithms, Gago- Vargas et al. (2005) concluded 
that the method of Briancon-Maisonobc is the best for the computation of s-parametric 
annihilators. In Levandovskyy and Martin-Morales (2008) we gave experimental results 
for the case / = /i and showed, that the algorithm of Briangon-Maisonobe is faster than 
the LOT method, which in turn is faster than the algorithm of Oaku and Takayama. 

Because of the structure of annihilators in the situation / = /i • . . . • /p, p > 1, basically 
the same principles stand behind the corresponding algorithms. 

Let s = (si, . . . , Sp), then a Bernstein- Sato ideal in K[s], which is defined as 



can be computed with the help of Ann/jj^,,] /* C D[s\. See Bahloul (2001) for algorithms. 
In contrary to the case / = /i, the ideal B{f ) need not be principal in general. However, 
it is an open question to give a criterion for the principality of B{f). Armed with such 
a criterion, one can apply the method of Principal Intersection 4.11 and thus replace 
expensive elimination above by the computation of a minimal polynomial. As in the case 
/ = /i it is an open question, which strategy and which orderings should one use in the 
computation of the annihilator in order to achieve better performance. 

3.5. Implementation 

Due to the comparison above, we decided to implement only Briangon-Maisonobe 

method for the (si, . . . , ,Sp)-paramctric annihilator Ann^ijs] C D[s\ in the case oip > 1. 

The corresponding procedure of dmod.lib is called annfsBMI. It computes both an- 
nihilator and the Bernstein-Sato ideal. 

We reported in Levandovskyy and Martin-Morales (2008) on several computational 
challenges, which have been solved with the help of our implementation. 

We use the following acronyms in adressing functions in the implementation: OT for 
Oaku and Takayama, LOT for Levandovskyy's modification of Oaku and Takayama 
(Levandovskyy and Martfn-Morales, 2008) and BM for Briangon-Maisonobe. Moreover, 



B{f) = {Km,n[: 



..]/r-----/;'' + (/i-----/p»n]K[si,...,Sp], 
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it is possible to specify the desired Grobner basis engine (std or slimgb) via an optional 
argument. 

For the classical situation where f — fi, there are Sannf sOT, Sannf sLOT, Sannf sBM 

procedures implemented, each along the lines of the corresponding algorithm. Moreover, 
there is a procedure Sannf s(/), computing Annj^j^j C D[s] using a "minimal user 
knowledge" principle. 

Example 3.2. We demonstrate, how to compute the s-parametric annihilator with 
Sannf s. This procedure takes a polynomial in a commutative ring as its argument and 
returns back a Weyl algebra of the type ring together with an object of the type ideal 
called LD. 
LIB "dmod.lib"; 

ring r = 0,(x,y),dp; // set up cominutative ring 

poly f = x~3 + y~2 + x*y"2; // define polynomial 
def D = Sannf s(f); // call Sannf s 

setring D; LD; // activate ring D, print Ann(f"s) 

==> LD [1] =2*x*y*Dx-3*x~2*Dy-y~2*Dy+2*y*Dx 

==> LD [2] =2*x~2*Dx+2*x*y*Dy+2*x*Dx+3*y*Dy-6*x*s-6*s 

==> LD [3] =x"2*y*Dy+y"3*Dy-2*x"2*Dx-3*x*y*Dy-2*y"2*s+6*x*s 

Note, that LD is not a Grobner basis but a set of generators. Computing a Grobner basis 
is done by groebner(LD) ;. In this case groebner returns the generators above and 2 
new ones. 

4. 6-functions with respect to weights for an ideal 

Let ^ w S ]R>o and consider the V-filtration V — {Vm \ Tn G Z} on D with respect 
to w, where Vm is spanned by {x^d^ \ —wa + w/S < m} over K. That is, Xi and di get 
w(ughts —Wi and Wi respectively. Note, that with respect to such weights the relation 
diXi = Xidi + 1 is homogeneous of degree 0. It is known that the associated graded ring 
©mez is isomorphic to D, which allows us to identify them. 

From now on we assume, that / is an ideal such that D / 1 is a holonomic module. Since 
holonomic D-modules are cyclic (e. g. Coutinho (1995)), for each holonomic D-module 
M there exists an ideal Im such that M = D/Im as £)-modules. 

Definition 4.1. Let ^ w € ]R>o. For non-zero 

p= Ca/3a;"5^ G D 

a,/3eNJ 

we put m = maxa^^{—wa + \ Cap 7^ 0} G M and define the initial form of p with 
respect to the weight w as follows: 

in(-,„,«,)(p) := 5Z Ca0X°'d^. 

ck,/5gNq: —wa-\-wfii—m 

For the zero polynomial, we set in(_,u, ,„)(0) 0. Additionally, we call the graded ideal 
in(-w,w)(^) ■— K • {iT^{-w,w)ip) \ P ^ 1} ttiG initial ideal of I with respect to the weight w. 
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Definition 4.2. Let ^ w G R>o <md s :— '^^—iWiXidi. Then in(_„, n K[s] is a 
principal ideal in K[s] . Its monic generator b{s) is called the global b-function of I with 
respect to the weight w. 

Theorem 4.3. The global 5-function of / is nonzero. 

We will give yet another proof of this theorem (Saito et al., 2000) in Section 4.2. 

Note, that by setting the weight vector in an appropriate way, one can compute b- 
functions of holonomic Z?-modulcs D/I, which are usually referred as 6- function for 
restriction, integration, localization etc. These special 6-functions play an important role 
in the computation of the corresponding restriction, integration, localization modules, 
sec Oaku (1997c); Saito ct al. (2000). 

Following its definition, the computation of the global 6-function of / with respect to 
w can be done in two steps: 

1. Compute the initial ideal I' of / with respect to w. 

2. Compute the intersection of /' with the subalgebra K[s]. 

We will discuss both steps separately, starting with the initial ideal. 

Computing the initial ideal 

In order to compute the initial ideal, the method of weighted homogenization has been 
proposed in Noro (2002). A more general approach on homogenization of differential 

operators can be found in Castro- Jimenez and Narvaez-Macarro (1997). 

Let u,v G R>o- The associative K-algebra D^^J^s^ is a G-algebra in the variables 
xi, . . . ,Xn,di, . . . ,dn,h which commute pairwise except for djXi = Xidj + 6ijh^^~^^^ . 
-^(u called the n-th weighted homogenized Weyl algebra with weights u, v, i.e. Xi and 

di get weights Ui and Vi respectively. 

For p = ^ CaffX^d^ G D one defines the weighted homogenization of p as follows: 

This definition naturally extends to a set of polynomials. Here, deg(„ „)(p) denotes the 
weighted total degree of p with respect to weights u, v for x, d and weight 1 for h. 

For a monomial ordering -< in D, which is not necessarily a well-ordering, we define 
an associated homogenized global ordering ^^''^ in d''^^^'^-^ by setting h ^'^^^ Xi, h di 
for all i and, 

p q if deg(„_^)(p) < deg(„^„)(g) 

or deg(„_^)(p) = deg(„_„)(g) and pi^^, ~< qi^^,. 

Note that for u — v = (1, . . . , 1) this is exactly the standard homogenization as in 
Saito et al. (2000). Analogue statements of the following two theorems can be found 
in Saito et al. (2000) and Noro (2002) respectively. Due to our different conception of 
Grobner bases (we require well-orderings), we give new proofs for them. 

Theorem 4.4. Let F be a finite subset of D and -< a global ordering. If G^'^^ is a Grobner 
basis of (iJ(„_„)(F)) with respect to -<^^\ then G^'^^i^^^ is a Grobner basis of {F) with 
respect to -<. 
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Proof. For any / G (F) with lm(iJ(„,^) (/)) = h^x°'d^, there exists g^''^ e G^''^ with 
lm(g('*)) = /I'^x'^a'^ satisfying lm{g'^^'>) \ lm(/). Then lm(5(''))|^^^ = x'^d^ \ x^d^ = 
lm(/)|,^^j, which proves the claim. □ 

Theorem 4.5. Let ^ be a global monomial ordering on D and -<[-w,w) the non-global 
ordering defined by 

x'^d^ -<[-w,w) if —wa + wjS < —w^ + wS 

or —wa + wjS = —wj + wS and x"d^ -< x'^d^. 

If G^''^ is a Grobner basis of H(^y^^y){I) with respect to ^l^"*^ then {in(_^^^^o)(5) I 9 S 
is a Grobner basis of 'm(^-w,wfi){H(u,v){I)) with respect to -<^"''. 

Proof. Let /' e in(_^^^^o)(-f^(ii,i;)(-^)) be (— 0)-homogeneous. There exist elements 
/ G H^uM^)'9 e GC') such that /' = in(_^ i„ o)(/) and Im (h) (g) \ Im (h) (/). 
Since /, g are {u, w)-homogeneous, we have 

Im (p) = lm^(M(in(_^ o)(p)) for p G {/.ff}, 

which finishes the proof. □ 

Summarizing the results from this section, we obtain the following algorithm to com- 
pute the initial ideal. 

Algorithm 4.6 (Initialldeal). 

Input: I C D such that D/I is holonomic, -< a global ordering on £), / w G K>0) 

M, V G 

Output: A Grobner basis G of in(_^^)(/) with respect to -< 
~^f^w w)'~ homogenized ordering as in Theorem 4.5 
G^'') := a Grobner basis of with respect to ^[^''^ 

return G = in(_„,.u,)(G('')|^^J 

4.2. Intersecting an ideal with a principal subalgebra 

We will now consider a much more general setting than needed to compute the global 
6-function. Let A be an associative K-algebra. We are interested in computing the inter- 
section of a left ideal J C A with the subalgebra K[s] of A where s G A \ K. We would 
like to find the monic polynomial b & A such that 

(6) = JnlK[s]. 

For this section, we will assume that there is a monomial ordering on A such J has a 
finite left Grobner basis G. 

Then we can distinguish between the following four situations: 

1. No leading monomials of elements in G divide the leading monomial of any power 
of s. 

2. There is an element in G whose leading monomial divides the leading monomial of 
some power of s. In this situation, we have the following sub-situations. 
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2.1. J ■ s C J and dimK(End^(A/J)) < cxd. 

2.2. One of the two conditions in 2.1. does not hold. 

2.2.1. The intersection is zero. 

2.2.2. The intersection is not zero. 
We now consider the first case. 

Lemma 4.7. If there exists no g € G such that Im(g') divides lm(s'^) for some k G Nq, 
then JnK[s] = {0}. 

Proof. Let 7^ 5 € J n K[s]. Then lni(5) = lm(.s'') for some fc G Nq. Since b e J, there 
exists g ^ G such that \in{g) \ lm(6) = lm(s'^). □ 

In the second situation however, we cannot in general state whether the intersection 
is trivial or not as the following example illustrates. 

Remark 4.8. The converse of the previous lemma does not hold. For instance, consider 
K[x, y] and J = {y^ + x). Then J n K[y] = {0} while {y^ + a;} is a Grobner basis of J for 
any ordering. 

In situation 2.1. though, the intersection is not zero as the following lemma shows, 
inspired by the sketch of the proof of Theorem 4.3 in Saito et al. (2000). 

Lemma 4.9. Let J • s C J and dimK(EndA(A/ J)) < oo. Then J n K[s] ^ {0}. 

Proof. Consider the right multiplication with s as a map A/ J — >■ A/ J which is a well- 
defined A-modulc cndomorphism of A/ .J as a — a' G J implies that (a — a')s G J • s C J, 
which holds by assumption for all a, a' G A. Since End^(A/J) is finite dimensional, 
linear algebra guarantees that this endomorphism has a well-defined non-zero minimal 
polynomial /i. Moreover, /i is precisely the monic generator of J H K[s] as fj,{s) = [0] in 
A/ J, hence /^(s) G J Pi IK[s], and deg(/Lt) is minimal by definition. □ 

Remark 4.10. In particular, the lemma holds if A/ J itself is a finite dimensional A- 
module. In the case where A is a Weyl algebra and A/ J is a holonomic module, we know 
that dimK(EndA(A/J)) < oo holds (e. g. Saito et al. (2000)). 

By the proof of the lemma, we have reduced our problem of intersecting an ideal with 
a subalgebra generated by one element to a problem from linear algebra, namely to the 
one of finding the minimal polynomial of an endomorphism. 

Proof of Theorem 4-3. Let w € R"q,/ C D such that D/I is a holonomic mod- 
ulo, J := in(_^u,^)(/) and s := '^^^I'WiXidi. Without loss of generality let 7^ p = 
13 (^a,ijx"d^ G J be {—w, i(;)-homogeneous. Then we obtain for every monomial in p 
by using the Leibniz rule 

x^'d^Xidi = x"'+^'d'^+^' + piX^d^ 

= {dix^'+^ - {at + l)xf')^d^ + pix'^d^ 

= {d,Xi - (a, + 1) + Pi)x"d^ = {x,d, -a, + Pi}x"d^. 
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Put m = —wa + wf3 for some term Ca.px^d^ in p where Ca.fi is non-zero. Since p is 
{—w, w)- homogeneous, m does not depend on the choice of this term. Hence, 

n n 

p ■ s = p ^ WiXidi = ^ '^{xidi -ai+ Pi)ca,i3X°'d'^ 

i=l i=l c«,/3 

n 

= s • p + ^ ^ Wi{-ai + Pi)ca,i3X°'d^ = {s + m) - pe J. 

i=l a,l3 

Therefore, J ■ s <Z J holds. Since D/J is holonomic (Saito et aL, 2000), Remark 4.10 and 
Lemma 4.9 yield the claim. □ 

If one knows in advance that the intersection is not zero, the following algorithm can 
be used for computing. 

Algorithm 4.11 (Principallntersect). 

Input: s e A, J c A a left ideal such that J n K[s\ ^ {0}. 

Output: b e K[s] monic such that J fl K[s] = (6) 

G :— a. finite left Grobner basis of J (assume it exists) 

i := 1 

loop 

if there exist ao,.. . , a,-! G K such that 

NF(s% G) + Y!j~Jo aj NF(sJ, G)=0 then 

return 6 := s' + X^^Cg '^i'^"' 
else 

i:=i + l 
end if 
end loop 

Note that because NF(s\ G) + X*^q NF(.s^ , G) = is equivalent to s*-|-X*~p ajs^ G 
J, the algorithm searches for a monic polynomial in K[s\ that also lies in J. This is done 
by going degree by degree through the powers of s until there is a linear dependency. 
This approach also ensures the minimality of the degree of the output. The algorithm 
terminates if and only if J n K[s] ^ {0}. Note that this approach works over any field. 

The check whether there is a linear dependency over K between the computed normal 
forms of the powers of s can be done by means of linear algebra. 

4-2.1. Applications 

Apart from computing global 6-functions, there are various other applications of Al- 
gorithm 4.11. 

4.2.1.1. Solving zero- dimensional systems Recall that an ideal I C K[a;i, . . . ,a;„] is 
called zero-dimensional if Klxi, . . . , Xn]/I is finite dimensional as a K- vector space. It is 
known (e. g. by Lemma 4.9) that in this case there exist ^ fi <E I D K[xi] for each 
1 < « < n, which implies that the cardinality of the zero-set of / is finite. 

In order to compute this zero-set, one can use the classical triangularization algorithms. 
These algorithms require to compute a Grobner basis with respect to some elimination 
ordering (like lexicographic one), which might be very hard. 
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By Algorithm 4.11, a generator of /nK[.r,;] can be computed without these expensive 
orderings. Instead, any ordering, hence a better suited one, may be freely chosen. 

A similar approach is used in the celebrated FGLM algorithm (Faugere et al., 1993). 
See also Noro and Yokoyama (1999) for a different approach. 

4-2.1.2. Computing central characters Let A be an associative K-algebra. The intersec- 
tion of a left ideal with the center of A, which is isomorphic to a commutative polynomial 
ring, is important for many algorithms, among other for the computation of the central 
character decomposition of a finitely presented module (cf. Levandovskyy (2005b)). In 
the situation, where the center of A is generated by one element (which is not seldom), 
we can apply Algorithm 4.11 to compute the intersection (known to be often quite non- 
trivial) without engaging much more expensive Grobner basis computations, which use 
elimination. 

Example 4.12. Consider the universal enveloping algebra of the Lie algebra SI2, A = 
U{sl2, K) = K(e, /, h I [e, /] = h, [h, e] = 2e, [h, f] = -2/). It is known, that over a field 
of characteristic 0, the center of A is K[4e/ + h"^ — 2h]. Consider a left ideal L and a 
two-sided ideal T, both generated by G = {e^^, /^^, - lOh^ + 9h} C A. Then consider 
A- modules Ml = A/L and Mt = A/T, which turn out to be finite-dimensional over K. 
We are interested in intersecting L.T with Z{A) and factorizing the output polynomial 
in one variable. The implementation of the Algorithm 4.11 in the library bfun.lib is 
described in Section 7. 

LIB "ncalg.lib"; LIB "central . lib" ; LIB "bfun.lib"; 

def A = makeUsl(2) ; setring A; // U(sl_2,Q) 

ideal Z = center (2); // generators of deg <= 2 

poly z = Z[l]; // we know there is just 1 generator 

ideal I = e-ll,f-12, (h-3)*(h-l)*h*(h+l)*(h+3) ; 

ideal L = std(I); // left GB of I 

vdim(L) ; // K-dimension of A/I 



// pretty-print factorization of polynomials: 
print (matrix(f actorize(vec2poly(imap(A, vT) ) , 1) ) ) ; 
==> z-3,z,z-15 

print (matrixCf actorize(vec2poly(imap(A, vL) ) , 1) ) ) ; 
==> z-3,z,z-440,z-8,z-48,z-168,z-15,z-99,z-120, 

z-255 , z-483 , z-575 , z+1 , z-399 , z-143 , z-195 , z-63 , 

z-80 , z-288 , z-360 , z-224 , z-323 , z-35 , z-24 

Note, that all the computations, thanks to Algorithm 4.11, were completed in a couple 
of seconds, while the Grobner-driven approach was still running after 20 minutes. 



==> 559 

vector vL = pIntersect(z,L) ; 
ideal T = twostd(I) ; 
vdim(T) ; 
==> 21 

vector vT = pIntersect(z,T) ; 
ring r = 0,z,dp; 



// L \cap K[z] 

// twosided GB of I 

// K-dimension of A/T 



// T \cap K[z] 

// commutative univariate ring 



13 



4-3. Intersecting an ideal with a multivariate subalgebra 

We now consider the case where we intersect J with the subalgebra K[s] = K[si, . . . ,Sr] 
of an associative K-algebra A for nonconstant, pairwise commuting Si, . . . , Sr & A. 

The following result is a consequence of a well-known characterization of zero-dimen- 
sional ideals. 

Lemma 4.13. The ideal JnK[s] is zero-dimensional if and only if for all 1 < z < r there 
exist /i S J such that lm(/i) = s^* for some di S Nq. 

Lemma 4.14. For a finite left Grobner basis G of J, 

GK.dim(K[s]) > GK. dim(]K[s]/(J n K[s])) 

> GK. dim(]K[s]/(i(G) n K[s])). 

Proof. For all / e J fl K[s] there exists g £ G such that lm{g) \ lm(/), which implies 
lm(5) G K[s] and thus, the claim follows. □ 

Note that the first incquahty is strict if and only if J n K[s] ^ {0}. 
We give a generalization of Algorithm 4.11 to compute a partial Grobner basis of 
J n ]K[s] up to a specified bound k gN. 

Algorithm 4.15 (intersectUpTo). 

Input: si,. . . ,Sr G A pairwise commuting, J c A a, left ideal, fc G N an upper degree 

boimd 

Output: a GB for J n K[si, . . . ,Sr] up to degree k 
G := a partial left Grobner basis of J consisting of elements up to degree k 

d:=0 
B:=% 

while d <k do 

Md := {s" I |a| < d} 

if there exist am G K, not all 0, such that ^ am NF(m, G) = then 

if J2 o.m'm ^ (B) then 

B:=BU{ "-mm} 

end if 
end if 

d:=d+l 
end while 
return B 

A couple of improvements can be made to speed up the computation time. 

If p G B with lm(p) = m has been found, any monomial which is a multiple of m can 
be discarded in the following iterations. 

Let G be a Grobner basis of J with respect to some fixed ordering -<. By using 
p G JnK[s] if and only if lm(p) G L(G)nK[s], one may disregard {m G | max^(m' G 
L{G)nMd) ^ m}. 

Further note that NF(m, G) = m, if m ^ L(G) n K[s]. 
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Using these improvements and ehoosing ^ to be a degree ordering and the elements in 
B to be monic, the output of the algorithm equals the reduced Grobner basis of J fl K[s] 
with respect to -< up to degree k. However, in general no termination criterion is known 
to us yet, that is apriori we do not know when we already have the complete needed basis 
of the intersection. Nevertheless, the termination is predictable if J fl K[s] is a principal 
ideal in K[s]. This situation often arises in the computation of Bernstein-Sato ideals, see 
Section 3.4. Moreover, another possibility for the algorithm to stop will be when the set 
of monomials we consider becomes empty on some step, which is the case if and only if 
J n ]K[s] is zero-dimensional. 

As one can see, the results above can be generalized by replacing the commutativity 
condition for a subalgebra S with the condition, that S* is a G-algebra in a K-algebra A. 
This and further generalizations will be studied in the next articles. Note, that under some 
extra requirements the algorithm will terminate after finally many steps without setting 
an explicit degree bound. Hence, in such cases a generally complicated elimination with 
Grobner bases can be replaced by much easier and predictable Grobner-free approach. 
The latter will, of course, allow to solve harder computational problems. 

As it was noted in Levandovskyy (2006), even the existence of a certain elimination 
ordering in G-algebras is not guaranteed. Consider the algebra B = K(a;, y \yx = xy+y"^). 
Then the ordering condition of Def. 2.1 says x > y must hold for any ordering. Hence, 
we cannot use Grobner basis in this G-algebra for computing the intersection of an ideal 
/ with the subalgebra K[a;], since the latter requires the use of ordering with x < y. One 
possibility would be to consider B as a, K-algebra with the ordering x < y modulo the 
two-sided ideal, generated by y^ —yx + xy. But this ideal has infinite two-sided Grobner 
basis, hence doing the elimination via passing to K-algebra setting is problematic, since 
it depends on the input ideal I. 

Despite these complications, it is obvious, that the preimage of an ideal in a subalgebra 
does exist. Hence, Algorithm 4.15 is indeed the only computational possibility to get some 
information about such a preimage. 

5. Bernstein- Sato Polynomial of / 

Let / e K[a;i , . . . , Xn] ■ One possibility to define the Bernstein-Sato polynomial of / is 
to apply the global 6-function for specific weights. 

Definition 5.1. Let B{s) denote the global ^-function of the univariate Malgrangc ideal 
If of / (cf. Section 3) with respect to the weight vector w = (1,0, ... ,0) € M""*"^, that 
is the weight of dt is 1. Then b{s) := B{—s — 1) is called the global b-function or the 
Bernstein-Sato polynomial of /. 

By Theorem 4.3, b{s) ^ holds. Moreover, it is well known that —1 is always a root 
of the Bernstein-Sato polynomial for nonconstant / and Kashiwara proved that all its 
roots are negative rational numbers (Kashiwara, 1976/77). 

The following version of Bernstein's theorem (Bernstein, 1971) gives us another option 
to define the Bernstein-Sato polynomial. 

Theorem 5.2 (Saito et al. (2000)). The Bernstein-Sato polynomial b(s) of / is the 
unique monic polynomial of minimal degree in K[s] satisfying the identity 

P • Z^+i = b{s) ■ for some operator P e D{s]. 
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Since P ■ f — b{s) e Ann(/^) holds, the Bernstein-Sato polynomial is the monic poly- 
nomial of minimal degree in K[s] that also lies in Ann(/*) -|- (/), hence b{s) is the monic 
generator of this intersection. 

Remark 5.3. There are the following choices for computing the Bernstein-Sato polyno- 
mial: 

(1) Compute a Grobner basis either of 

(a) J = in(_u,,u,)(//), which amounts to 1 Grobner basis computation in £>„ or 

(b) J = Ann(/*) -I- (/), which requires 2 Grobner basis computations in D„[s]. 

(2) Intersect J with K[^] where ^ is either s or WiXidi 

(a) the classical elimination-driven approach (needs 1 tough Grobner basis compu- 
tation) or 

(b) using Algorithm 4.11 with no Grobner basis computation. 

It is very interesting to investigate the new approach for the computation of Bernstein- 
Sato polynomials, arising as the combination of the two methods 

(1) Ann(/*) via Briangon-Maisonobe (cf. Section 3 and Levandovskyy and Martin- 
Morales (2008)), 

(2) (Ann(/'*) + (/)) n K[s] via Algorithm 4.11. 

For the computation of in(_„, „,)(//) using the method of weighted homogenization as 
described in Section 4.1, the following choice of weights is proposed in Noro (2002) for 
an efficient Grobner basis computation: 

u = (degji(/),Mi, . . .,Un), 

V = (l,degji(/) - ill + 1,. . .,degs(/) -u„ + 1) 

such that the weight of t is deg^(/) and the weight of dt is 1. Here, u G M"g is an 
arbitrary vector and deg^(/) denotes the weighted total degree of / with respect to u. 
The vector u may be choosen heuristically in accordance to the shape of / or by default, 
one can set u = (1, . . . , 1). 

6. Enhancements to steps of algorithms 
6.1. Enhanced computation o/ Ann£)[g](/*) 

Consider the set of generators G := {fdt + s, {fidt + di \ 1 < i < n} of an ideal 
J, coming from the Briancon-Maisonobe method. According to the latter, we have to 
eliminate dt from J, that is to compute Jfl D„[s] = Ann£)[s] (/**). 

Since any element h from J has a presentation as 

n n n 

h = ao{fdt + s) + ^ai{fidt + di) = {aof + ^aifi)dt + {aos + ^aidi), 

i=l i—1 i=l 

then for all (ao, ai, . . . , a„) € Syz({/, /i, . . . , /„}) n IK[.t, .s]"+-^ we obtain that oqs + 
Z]r=i "i^i ^ ^ ri Dn[s]. 

Moreover, it is known, that indeed the above elements generate the ]K[a;]-submodule 
of all the elements in J n £>n[s], which total degree in di does not exceed 1. 

Consider the set Tf = {/, /i, . . . , /„} C K[s] C Dn[s], a left ideal _D„[s]T/ C Dn[s] and 
an ideal K[x]T/ C K[x]. Then a Grobner basis of ]K[a;]T/ is a Grobner basis for I?„[s]T/ 
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as well. Denote by Sf a set of generators of the module Syz(T/) C K[a;]"+-'^. By e. g. 
generalized Schreyer's theorem (Levandovskyy, 2005a), it follows that the module of left 
syzygies LeftSyz^i^j^j (T/) = Dn[s]Sf. 

Let ^1 be a monomial module ordering on K[a;]""'""'^, which is a position-over-term 
ordering, which gives preference to the 1st component. Since degree of / is always by 1 
bigger than the degree of ^ , the cofactors to / have respectively smaller degree. 

Algorithm 6.1 (Sannf sBMSyz). 
Input: / e K[x] 
Output: Ann£,„[^](/^) 



Remark 6.2. One of major difficulties in the computation of Grobner basis (especially 

with respect to an elimination ordering) is the need to compute numerous intermediate 
polynomials (of usually high degree and with big coefficients) in order to come to a 
polynomial in the answer, which is often of small degree with coefficients of moderate 
size. Actually the set of generators Sa, which we compute in the syzygy-driven algorithm, 
generates already a part of the answer, though the corresponding ideal is, in general, not 
yet the complete answer. 

Computing a Grobner basis of Sa and adding it to the original set of generators G 
allows to avoid at first place the discovery of elements of Sa in the Grobner basis computa- 
tion of GuSa and hence allows to decrease the number of intermediate unpleasant polyno- 
mials, which are needed in such computation. This is important, since in the answer there 
are no polynomials of degree zero with respect to di, that is Annjj^[s]{f^) r\K[x,s] = 
due to the fact, that the only element from the ring ]K[a;, s], annihilating is zero. 
Hence with Sa we add the set of elements of smallest possible total degree in di that is 
of degree 1. Such elements are, in general, very hard to compute via the Grobner-driven 
elimination. 

However, it is very interesting to derive conditions, under which the above algorithm 
is more efficient than the one of Briangon-Maisonobe. We observe that it is not true for 
a couple of examples. See section 7.1. 

6.2. Enhanced computation of bf{s) 

In the following Lemma we collect folklore results and supply them with short proofs 
for the completeness of exposition. 

Lemma 6.3. Let, as before, / e K\x] \ {0}. 



(1) V 1 < i < n we have - s-§^ G Annuls] and -§§-dj - G AnnD[s] /' 

(2) For / G K, bf{s) = 1. For / G k[x] \ {K}, {s + 1) \ b}{s). 
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(3) 

(^> = (A„„r + tf,^,.^..|i»nKW 

(4) if K = K, then V{f) is smooth implies 6/(s) = s + 1. 

df . . . 

Proof. Wc use shortcut fi :— — ^. Consider surjoctive K-alg homomorphism tTq : D\s\ — > 

D, s H- ?■ a and apply it to the inclusion Dls]{P{s)f — i>{s)) C Annxjj-^] f^. Then we have 
an inclusion 

7r„((P(s)/ - 6(s))) = {P{a)f - 6(a)) C (Ann^[,] /«) |,=„C Ann,, /« 

(1) Direct calculation. 

(2) By using 7r_i from above, we obtain P(— 1) = &(— 1)/~^ modulo Ann£) /° = 
{di, . . . ,dn)- Hence -P(— 1) = p{x) € K[a;] and p(a;) = 6(— 1)/~^, which can be 
true only in two cases: 

1. / e K*, then P{s) = /"i G K and 6(s) = 1, 

2. / ^K*, then 6(-l) = 0. 

(3) Let us write P{s) = J2i Pidi+Po for Pq € K[x, s] and Pj € -D[s]. Computing modulo 
Ann£,[s] and using (1) we can present P(.s)/ = Pof + Y^i Pidif = Pn.f + {s + 
1) Ylii Pifi- By (2) 6(— 1) = 0, hence by specializing s to —1 in Bernstein's equation 
we get P(-l) • 1 = b{-l)f-^ = 0. Thus P(-l) e Annc[,](l) = (91,..., 5„). In 
particular, Po(— 1) = and hence s + 1 | Pq G K[x, s]. Moreover, 

E^^/^ + ^Z-S^eAnn^^r 

i 

and the claim follows. 

(4) Since / = I{Sing{V{f))) = (/, /i, ■ ■ ■ , /«) C lK[a;], smoothness takes place when 
1 e /, hence by (3) we have 1 G (j^) and thus, b{s) = s + 1. □ 

Lemma 6.4. For a fixed algorithm, which computes Ann^ijs] (/"), let us consider two 

ideals Ji = Ann,, [,](/*) + (/) and /a = Ann£,[s] (/^) + (/, . . . , ^) (note that /i, /a C 
£'[s]). There are two following algorithms, which take an ideal and a polynomial as input 
and return Bernstein-Sato polynomial, namely 

Algorithm 1. b{s) = pIntersect(/i , s). 

Algorithm 2. b{s) = (s + !)• p1ntersect(/2, s). 
Then the Algorithm 2 is more efficient than the Algorithm I. 

Proof. Performing the principal intersection, the Algorithm 2 will compute one normal 
form less (of an element of high degree) than the Algorithm 1. Moreover, the normal 
forms in Algorithm 2 are taken with respect to a bigger ideal, what makes respective 
computations easier as well. By Lemma 6.3 (1) we know that fdi — sfi G Ann£3|-^](/''). 
Hence these elements can be reduced to (s + 1) • fi in 7i = Ann£)[s] (/*) + (/). Meanwhile 
in I2 = Annc[s](/*) + (/, /i, . . . , /„) we reduce fdi — sfi automatically to zero. Note, 
that indeed (s + l)/2 C I\ C I2 holds and hence, in the process of computing a Grobner 
basis of I\ (Algorithm I), the operations with commutative elements of the kind (s + l)/i 
will in general keep the factor (s + I), thus operating with larger polynomials of higher 
degree. Hence the claim. □ 
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6.3. Enhanced computation of normal forms 

When computing normal forms of the form NF(s*, J) like in Algorithm 4.11 we can 
speed up the reduction process by making use of the previously computed normal forms. 

Lemma 6.5. Let ^ be a K-algebra, J C A a, left ideal and let f G: A. For i e N put 
n = NF(/% J), Qi = f"- - n e J and Ci = j^j^T^ provided nqi ^ 0. For nqt = we put 
Cj = 0. Then we have for all z £ N 

r,+i = NF(/r„ J) = NF([/^ - r„ n],. + r.n, J). 

Proof. It holds that /'+^ = fqi + fri fri, which shows the first equation. On the other 
hand, /*+^ = Qif + rif = qi{qi+ri)+ri{qi+ri) = qiqi+qiri+riqi+rtn qiri+rin = 
(/' — ri)ri + rjri [/* — ri,ri]cj + r-jri, which proves the second equation. □ 

As a direct consequence, we obtain the following result for some K-algebras of special 
importance. 

Corollary 6.6. If ^ is a G-algebra of Lie type (e. g. a Weyl algebra), then 

n+i = NFifn, J) = NF([r - n, n] + nn, J) holds. 

If A is commutative, we have 

n+i = NFinruJ) = m{r\+\j). 

Note, that computing Lie bracket [f,g] both in theory and in practice is easier and 
faster, than to compute [f,g] as f ■ g — g ■ f, see e. g. Levandovskyy and Schonemann 
(2003). 

Remark 6.7. We work on enhanced algorithms for the computation of the Bernstein 
operator from Theorem 5.2 as well and will report on the progress in forthcoming articles. 

7. Implementation 

In Noro (2002), M. Noro proposed methods of modular change of ordering and modular 

solving of linear equations to be used in his approach, which is based on a kind of 
Algorithm 4.11. In our implementation we decided to develop, test and enhance first 
purely characteristic methods, thus having the possibility to adjoin modular methods 
later. 

For the computation of 6-functions and Bernstein-Sato polynomials, we offer the fol- 
lowing procedures in the Singular library bfun.lib: 

bfct computes in(_^ „,)(//) using weighted homogcnization with weights u,v for an 
optional weight vector u (by default u = (1, . . . , 1)) as described above, and then uses 
Algorithm 4.11 with the enhancement from Corollary 6.6 for the intersection, where the 
occuring systems of linear equations are solved by means of linear algebra. 

bf ctSyz computes in(_^^^)(//) as in bfct and then uses Algorithm 4.11, where the 
linear equations are treated as polynomial ones and then solved by computing syzygies. 

bf ctAnn computes Ann(/®) via Algorithm 6.1 and then computes the intersection of 
Ann(/*) + (/, . . . , ^) with K[s] analogously to bfct. 
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bf ctOneGB computes the initial ideal and the intersection at once using a homogenized 
elimination ordering (see also Hartillo-Hermoso (2001)). 

For the global 6-function of an ideal /, bf ctldeal computes in(_u,, ,„)(/) using standard 
homogenization, i. e. weighted homogenization where all weights are equal to 1, and then 
proceeds the same way as bf ct. Recall that D/I must be holonomic as in Saito et al. 
(2000). 

All these procedures work as the following example illustrates for bf ct and the hy- 
perplane arrangement xyz{y — z){y + z). 
LIB "bfun.lib"; 
ring r = 0, (x,y,z) ,dp; 
poly f = x*y*z*(y-z)*(y+z) ; 
bfct(f) ; 
==> [1]: 
==> _[!]=-! 
==> _[2]=-5/4 
==> _ [3] =-3/4 
==> _ [4] =-3/2 
==> _[5]=-l/2 
==> [2]: 

==> 3,1,1,1,1 
7.1. Comparison 

We use the polynomials in Table 1 for test examples, where we measure the total 
running time of each call to a system in a batch mode. In this time the initialization of a 
system, loading of an example file, the actual computation and the writing of an output 
are included. 

Table 1. Examples 



Example 


Input 


ab23 


{z^ + w^){2xz + 3yw^) 


cnu6 


{xz + y){x^ -y^) 


cnu7 


{xz + y)(x'^ - j/^) 


tt43 


x^ + y* + z*{xyzf 


xyzcusp45 


{xz + y){x* + y^) 


uwl8 


xyz{x ~ z){z — y){y + z){2x + 2y — z) 


uw22 


xyz{x + z){x- y){y - z){y + z) 


uw27 


xyz{x + y){-x + 2y + z){x + y + z){y + z) 


uw28 


xyz{x + z){-x + y + z)(x + y){y + z) 


uw29 


xyz{x + z){x + j/)(-3.x- - J/ + 2z){y + z) 


uw30 


xyz(x - 2z){-x + y + z){x - y){y + z) 



20 



The running times in the; tables below are given in "[hours[h]:]niinutc!s:seconds" format. 
We use the shortcuts when we have stopped the process after the time t and when 
the process ran out of memory after the time t. 

The tests were performed on a machine with 4 Dual Core AMD Opteron 64 Processor 
8220 (2800 MHz) (only one processor could be used at a time) equipped with 32 GB 
RAM (at most 16 GB were allowed to us) running openSUSE 11 Linux. 

Wc first request the computation of Ann/j^jsj (/*) and the Bernstein-Sato polynomial 
comparing the different algorithms from Section 6. We use the notation from Lemma 6.4. 

Table 2. Comparison of the algorithms from Section 6 





Anno^fs 


i(r) 


Bernstein-Sato polynomial 








Sannf sBMSyz based 


Sannf sBM based 


Example 


Sannf sBMSyz 


Sannf sBM 


Alg. 1 


Alg. 2 


Alg. 1 


Alg. 2 


ab23 


0:01 


0:02 


0:05 


0:03 


0:06 


0:03 


cnu6 


0:01 


0:01 


0:01 


0:01 


0:01 


0:01 


cnu7 


0:09 


0:12 


0:11 


0:11 


0:19 


0:12 


tt43 


0:01 


0:01 


0:03 


0:01 


0:03 


0:01 


xyzcusp45 


0:56 


1:22 


1:30 


1:16 


1:20 


1:10 


uwl8 


3:38 


0:04 


17:36 


12:38 


16:24 


11:15 


uw22 


2h:04:01+ 


0:04 


2h:17:36t 


2h:16:55+ 


2h:28:07''' 


2h:28:53+ 



Further, we compare our implementations for the computation of the Bernstein-Sato 
polynomial with the existing ones in the computer algebra systems Risa/Asir and 

MACAULAY2. 

We have used Risa/Asir version 20071022, Macaulay2 version 1.1 with version 1.0 
of Dmodules .m2 and SINGULAR 3-1-0 with bfun.lib version 1.13. 

We would like to stress, that in our implementation of bfun.lib we have restricted 
ourselves to the use of characteristic zero methods, in order to see what can we achieve 
with them. The implementation of AsiR by M. Noro (Noro, 2002) uses the methods in 
prime characteristic, which can be applied to our implementation as well. However, the 
values in the table above indicate, that the difference in timings is not devastating for 
our cause. 

As the timings in Table 3 suggest, the approach via the initial ideal seems to be 
specially well suited for hyperplane arrangements, while it looks like that the performance 
of the annihilator based method is better for other kind of input (we took non-quasiho- 
mogeneous singularities). See Walther (2005) for details about generic arrangements. 
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Table 3. Comparison of different systems 







Asm 


MACAULAY2 


Singular 


Example 


bfimction bfct 


globalBFimction 


bfct 


bfct Ann 


ab23 


0:23 


0:17 


0:27 


0:17 


0:04 


cnu6 


1:39 


0:54 


14:03 


0:01 


0:01 


cnu7 


7:32 


4:46 


4h:03:39^ 


0:06 


0:20 


tt43 


0:07 


0:05 


0:05 


0:17 


0:01 


xyzcusp45 


1:52 


1:10 


4h:18:35 


3:05 


3:01 


uwl8 


7:22 


29h:35:54^ 


4h:08:16'' 


6:21 


12:27 


uw22 


2:12 


4h:04:05'' 


4h:01:43^ 


2:24 


2h:42:02 


uw27 


2:37 


3h:05:14^ 


llh:45:18'' 


4:40 


6h:55:35^ 


uw28 


1:36 


10h:23:40^ 


3h:03:00^ 


3:10 


3h:03:32>< 


uw29 


1:48 


3h:51:14^ 


10h:23:42'' 


2:52 


3h:01:30'' 


uw30 


1:58 


5h: 14:18'' 


3h:06:57^ 


3:09 


3h:00:13'< 



8. Bernstein-Sato Polynomial for a Variety 

In the paper of Budur, Mustafa and Saito (Budur et al., 2006), using the theory of In- 
filtrations of Kashiwara (Kashiwara, 1983) and Malgrange (Malgrange, 1983), the theory 
of the Bernstein-Sato polynomial of an arbitrary variety has been developed. We present 
here the construction of the Bernstein-Sato polynomial of an afHne algebraic variety. 

Given two positive integers n and r, for the rest of this section we fix the indices 
i,j,k,l ranging between 1 and r and an index m ranging between 1 and n. 

Let / = (/i, . . . , /r) be an r-tuple in K [a;] '".Consider a free K[x, s, j]-module of rank 
one generated by the formal symbol and denote it by M = K[x, s, j] • Here, s = 
(si, . . . , Sr), J = \j and = /^^ ■ • ■ ff.''. Moreover, we denote by K{S) the universal 
enveloping algebra U{qI^), generated by the set of variables S = (Sy), i,j = l,...,r 
subject to relations: 

[sij,Ski] = SjkSii - SiiSkj- 

Then, we denote by Dn{S) := £>„ (8>k K(S'), which is a G-algebra of Lie type by e. g. 

Levandovskyy and Schoncmann (2003). 

The module M has a natural structure of left Z)^ (<S')-module when the variables 
act in the following way {i < j): 

Sij»{G{s)-n = Si-G{s + ej-ei)^-f e M, 

Ji 

where G{s) is an element in K[x, s, j] and ej stands for the j-th basis vector. 

One can easily observe that the action of Sa on M coincides with the multiplication 
by Si from the left. 
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Following the ideas by Malgrange, one can also consider M as a Dn{R) 0k Dj-iT)- 
module, with T = Wyt],t = (ti, . . . , i^), dt = {dti, . . . , dtr) and the action 



dti • (G(s) • n = -SiG{s -ei)j:- f\ 
Observe that the action of above corresponds to the action of —dti ■ tj- 



(1) 



Theorem 8.1 (Budur, Mustata, and Saito (2006)). For every r-tuplc / = (/i. .... /, ) G 
M.[xY there exists a non-zero polynomial in one variable b{s) G K[s] and r differential 
operators Pi (5), ... , PriS) G Dn{S) such that 

r 

Y,Pk{s)fk-r = b{s, + --- + sr)-r e M. (2) 

fe=l 

The Bernstein-Sato polynomial bf{s) of / = (/i, • • • , fr) is defined to be the monic 
polynomial of the lowest degree in the variable s satisfying the equation (2). ft is demon- 
strated in Budur et al. (2006), that every root of the Bernstein-Sato polynomial is ra- 
tional. Let I be the ideal generated by /i , . . . , /j. and Z the (not necessarily reduced) 
algebraic variety associated with / in K". Then it can be verified that bf{s) is independent 
of the choice of a system of generators of /, and moreover that bz{s) = bf{s — codimZ+l) 
depends only on Z. For instance, the Bernstein-Sato polynomial of f{x, y) G K[x, y] and 
Bernstein-Sato polynomial of the variety defined by the ideal {f(x,y),z) coincide. How- 
ever, due to the codimension, there is a shift between 6/(s) and ^'(/(a;,j/),2)(s)- 

Now, let us denote by Annjj^(^s){f^) the left ideal of all elements P{iS) G Dn{S) such 
that P{S) • = 0. We call this ideal the annihilator of in Dn{S). Prom the definition 
of the Bernstein- Sato polynomial it is clear that 

(AnnB„(s)(r ) + (/i, • • • , /r)) n K[si + . • • + s,] = {bf{si + . . . + s,)). 

Since the final intersection can be computed with the Principal Intersection method 
4.11, the above formula provides an algorithm for computing the Bernstein- Sato polyno- 
mial of afhne algebraic varieties, once we know a Grobner basis of the annihilator of /* 
in Dn (S) . The rest of this section is dedicated to the solving of this problem. 

8.1. The annihilator of f "^ in D„(5) 

Consider the generalization of Malgrange's ideal If associated with / = (/i, . . . , fr), 

I ^ dfi l<i<r 

If = {ti-fi, dm + ^g^dt^ - - 



C Dn(t,dt) 

1 < m < n 



Here we give a computer-algebraic proof to the following Lemma, whose assertion is 
expected as in Saito et al. (2000) (for instance). 

Lemma 8.2. If is a maximal ideal in Dn{R) ®k Dr{T) and // = ATxn.jj^(^ii-^(^^D^(j-) 
Proof, {ti - fi) • f = fif - fif = 0. For any m, observe that 

dm • (/r •...■/-)= ^ 5„ . (/f )(/- ■ . . . • fnif-T' = E ^^^^(/^ • • • • • /") 
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Since dtk acts on by the multiplication with — Sj/j ^, the generators of the second 
type annihilate so If C AmiD^(R)®^Dr(T) 

By Lemma 3.1, the set of generators of // is the same as the set S\ in the Lemma 
and hence there is a monomial ordering, such that Si is a Grobnor basis. In particular, 
If is a proper ideal. The set of leading monomials of 5*1 is then L = {tj,dm}- Since any 
monomial ordering on can be presented as weighted degree ordering with the 

weight vector w with strictly positive entries (see e. g. Bueso et al. (2003)), we see that 

GK.dim(D„(i?) Ok Dr{T))/If = GK.dimK[{tj,dtj,Xm,dm}]/{L) = r + n. 

Assume the left ideal // is not maximal, then there exists p^If, such that // C 
If + (p) c Dn{R) Ok Dr{T). In particular, \m{p) docs not include the elements of L 
above. If If + (p) is a proper ideal, its set of leading monomials strictly includes L and 
has at least one element more. But then the dimension argument as above shows, that 
GK. dim(Z)„(i?) Ok Dr{T))/{If + (p)) < r + n, what contradicts Bernstein's inequality. 
Hence // is maximal and it is equal to the annihilator. □ 

Theorem 8.3. Let / = (/i, . . . , /,.) be an r-tuplc in K[xY and Dn{dt, S) the K-algebra 
generated by dt and S with the corresponding non-commutative relations. Then the 
following ideal of Dn{S) coincides with the annihilator of /* in Dn{S): 



Dn{dt, S) (^Sij + dtifj , dm+Y^ ^^^^^ 

k—l 



1 < i, j < r 
I < m < n 



nDr,{S). 



Proof. Let (j> : Dn{S) ^ Dn®¥^Dr{T) be the K-algebra homomorphism given by 4){sij) = 
—tjdti — 6ij and 0(P) = P for all P in D„. In view of Lemma 8.2, we observe that 

Annn^^s){n = Dn{S) nlf= 4>-\lt)- 

The morphism (f) can be written as (/> = ® 99, where : IK(5) ^ Dr{T) = K{t, dt), 
Sij I— >■ —tjdti — Sij. Thus we can apply Proposition 2.3 to obtain that Ann^j^^s) (/*) = 
{E'l^ + E'lf) n £»„(6'), with I^ = {{sij + tjdti + Sij I 1 < i,j < r}) and E' = 
Dn{{tj,dtj,Sij}) subject to relations 

[sij,Ski] = SjkSil - SiiSkj, [dtk,tj] = djk, [sij,tk\ = -Siktj, [sij,dtk] = Sjkdti. 

By Theorem 2.2, E' is a G-algebra, if there exists an elimination ordering for 
{dti, . . . , dtr} on Dn{dt, S), obeying the conditions 

\m{SjkSii - SiiSkj) < SijSki, tj < Sijti, and dti < Sijtj. 

It is clear, that such orderings exist. 

Now, we proceed with the elimination of {ti,dti \ 1 < i < r} from (// + I,p) in E' . 
By taking a monomial ordering with the property {tj} » {a;,}, {di, Sij} » {xi, dtj}, we 
start with eliminating {tj} first. 

By Lemma 3.1, the generators Gi of // form a Grobner basis. The ideal lip in the 
current situation is generated by G2 = {sij + tjdti + Sij}. In order to prove, that Gi UG2 
is a a Grobner basis, we apply the generalized Product Criterion (Lemma 2.4). At first 
we apply reduction process by Gi, thus obtaining G'2 = {sij + fjdti}. Then 

[sij + fjdti, Ski + fidtk] = Sjk{sii + fidti) - Sii{skj + fjdtk) 
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which clearly reduces to zero. The next kind of pairs 

[sij + fjdti,tk - fk] = -5iktj + 5ikfj = -Sik{tj - fj) 
again reduces to zero. It remains to consider 

[sij + fjdti, dm + ^ -^-^^tk] = ^ jr^Sjkdti - dti[dm, fj] = 0, 
k=i "^"^ k=i ™ 

since [dm,fj] = 

Hence, GiU G'2 is a a Grobner basis and hence, by the Elimination Lemma, G3 = 
(Gi U G2) \ {tj — fj} is a Grobner basis of {If + 1^) fl Dn{dtk, Sij). Thus, it follows, that 

Annnis){fn = (sij + dUfj ,dm + J2 P^^*^) n □ 

fe=i "^"^ 

Indeed, the result we have proved is a natural generalization of the algorithm for com- 
puting the annihilator of in Dn[s] (cf. 3.2) given by Briancon-Maisonobe in Briancon 
and Maisonobe (2002) . Finally, the algorithm for the computation of Ann^i^ (^s) f looks 
as follows: 

Algorithm 8.4 (SannfsVar). 
Input: / = (/i, . • . , /r), an r-tuple in K[xY 
Output: {G\{S), . . . ,Ge{S)}, a Grobner basis of AnnD^(^s){f^) 
Let Dn{dt, S) be the algebra in Corollary 8.3, with non-commutative relations 

[di,Xi] = 1, [sij,dtk] = Sjkdti, [sij,Ski] = SjkSii - SuSkj- 

Ji := {{s,, + dUf, \ l<i,j<r}) 

J2 ■■= {{dm + ELi it^t, \l<m<n}) 

J:=Ji + J2 >JCDn{dt,S} 
H := G.B. of J w.r.t. a compatible elim. ordering for dti, . . . , dtr 
HnD^iS) =: {G,{S),...,G,{S)} 
return {Gi{S), . . . ,Ge{S)} 

8.2. Elimination orderings in Dn{dt,S) 

One of the bottlenecks of the presented algorithm for computing the Bernstein- Sato 
polynomial for varieties is to calculate the corresponding annihilator. An elimination 
term ordering for {dti, . . . , dtr} in Dn{dt, S), which has 2n + r + r^ variables, has to be 
considered. In addition, due to the structure of the G-algebra, this ordering < has to be 
chosen with the following extra restrictions 

dti < Sijdtj, lm{5jkSii - SuSkj) < SijSki 

for all indices i,j,k,l where the expression makes sense. The efficiency of the method 
strongly depends on the selected ordering. Therefore it is worth analysing it in detail. 



25 



Assume that < is such an ordering and let us consider the first two rows of the matrix 
representing the ordering in this way. 



dti ■ ■ ■ dtr S X dx 



Pi ■ 


• Pr 


a b 


c 




■ Ir 


a P 


7 


<' 



The vectors a, b and c must be zero, since < is an ehmination ordering for {dti}. The 
conditions dti < Sijdtj, imply pi < pj for all Thus all pi, . . . ,Pr are equal and can 
be taken as 1. 

From computational point of view, since the variables {s^} do not commute with 
{dti}, these two blocks must be together in the elimination ordering, namely /3 = 7 = 0, 
otherwise Grobner bases computation may be slow. 

In the implementation we have taken an = 2 and atj = 1 for i 7^ j, and q = 0. 
However, in some examples we have observed that lexicographical orderings are also 
useful, see Example 8.8 below. 

In this section wc have described an algorithm for computing the Bernstein- Sato poly- 
nomial of affine algebraic varieties without any homogenization but passing through the 
computation of the annihilator of /" in Dn{S). Now, other methods are ilustrated. 

8.3. Another approach 

As Budur et. al. point out in (Budur et al., 2006, p. 794), the Bernstein-Sato polyno- 
mial for varieties coincides, up to shift of variables, with the 6- function in (Saito et al., 
2000, p. 194), if the weight vector is chosen appropriately. Let us describe this algorithm 
more carefully. 

Let // = Ann£)^^(_a()(/'') be the Malgrange ideal associated with / = (/i, . . . , /r) and 
consider the weight vector w = ((0, . . . , 0), (1, . . . , 1) e Z" x which gives weight to 
dm and weight 1 to dti. Consider the F-filtration V = {Vk \ k € Z} on Dn{t,dt) with 
respect to w, where Vk is spanned by {t" ■ dt^ \ —\a\ + |/3| < k} over K. Note that the 
associated graded ring ®/cezVfe/Vfe_i is isomorphic again to the (n + r)-Weyl algebra 
Dn{t,dt) and the homogeneous parts are the following. 

'Dn{t^-dtj)dP, |/3| = fc>0; 
Vk/Vk-i = lD„{t,-dt,), fc = 0; 

Dn{ti ■ dtj)t°' , -|a| = fc<0. 

Denote by B{s) the 6-function of the holonomic ideal // with respect to w. Recall that 

B{s) is the monic generator of the ideal in(_^„ u))(//) n 'K\t\dt\ + h trdtr]- 

As in the classical case, i.e. r = 1, the following result holds. 

Lemma 8.5. 6/(s) = (-1)<i^s^WB(-s - r). 

Proof. Consider Pi{S), . . . , Pk{S) G Dn{S) differential operators satisfying the func- 
tional equation X]fc=i Pk{S).fk • = bf{si + ■ ■ ■ + s,.) • Then fe/(si + • • • + s,.) — 
J2k=i Pk{S)fk is an element in Ann^,^^^^ (/*) and hence applying the Mellin transform. 
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or equivalcntly making the substitution i-)- —tjdti — Sij, one obtains the following 
element in //. 

r 

bf{-tidti trdtr -r) - Pk{-tjdti - Sij)fk e // 

k=l 

Modulo // the polynomials fk in the above expresion can be replaced by tk, since tk — 
fk S //. Finally, taking initial parts one concludes that bf{—t\dt\ — ■■ ■ — trdtr — r) € 
iii(-tu.ti))(-^/)i which means that B{s) divides bf{—s — r). 

Conversely, by definition there exists a differential operator P{t,dt) £ If C Dn{t,dt) 

such that B{tidti H h trdtr) = in(_TO_^)(P(t, dt)). In particular P{t, dt) has F-degree 

zero and hence it can be decomposed into y-homogcncous parts as follows 

P{t, dt) = B{tldtl + ■■■+ trdtr) + Qcx{tidtj)t". 

\a\>l 

Since P{t,dt) G //, making left reduction of P{t,dt) with respect to {ti — fi} we arrive 
at B{tidti H h trdtr) + Y.\a\>iQoc{tidtj)f°' & If r\ Dn{tidtj). After applying the 

substitution tidtj —Sji — Sij, wc conclude that B{—si — ... — — r) belongs to the 
ideal Ann/3^^5) (/*) + (/i, . . . , fr) and the proof is complete. □ 

Algorithms for computing this 6-function, which use the homogenization technique in 
the Weyl algebra, are given in Section 4, see also Saito et al. (2000). We describe the 
complete algorithm for computing Bernstein-Sato polynomials using initial parts. 

Algorithm 8.6 (bf ctVar). 

Input: / = (/i, . . . , fr), an r-tuplc in K[a;]''; Z. variety associated with / 
Output: bz{s) = bf{s — codimZ + 1), Bernstcin-Sato polynomial of Z 
Let Dn{t, dt) = Dn <8) Dr be the (n + r)-Weyl algebra. 

I ■■= {{i^ - . + E;=1 > I ^ Dn{t,dt) 

w := ((0,...,0),(1,...,1)) G Z" X Z"^ 

J := InitialIdeal(J, w) > Algorithm 4.6 

s := —{dti • ti + . . . + dtr ■ tr) 

b{s) := plntersect(s, J) > Algorithm 4.11 

return b{s — codim Z + 1) 

As for the Bernstein-Sato polynomial of a polynomial (indicated in Remark 5.3), we 
have two different methods for the computation of Bernstein-Sato polynomial of an affine 
algebraic variety, namely 

• minimal polynomial for si -|- . . . -|- Sr in D/ in(_tu „,)(//), 

• (Ann + (/i, . . . , fr)) n IK[si + . . . + Sr], where the intersection can be done rather 
with the 4.11 method, than by using Grobner basis elimination. 

It is important to investigate the connection of these methods and especially their 
applicability to different classes of varieties. Our experience shows, that no method is 
clearly superior to the other one in general. Thus it is desired to have both of them in 
any package for D-modules. 

Remark 8.7. Very recently, the authors have realized another approach for computing 
Bcrnstcin-Sato polynomials for varieties. In Shibuta (2008), Shibuta modifies the def- 
inition of Budur-Musta^a-Saito's Bernstein- Sato polynomials to determine a system of 
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generators of the multiplier ideals of a given ideal. Then he obtains an algorithm for 
computing Bernstein-Sato polynomials, which gives an algorithm for computing multi- 
plier ideals and jumping coefficients. His methods are based on the theory of Grobner 
bases in Wcyl algebras and corresponds to the natural generalization given by Oaku and 
Takayama, hence they need homogenization techniques. 

We conclude this section showing several examples calculated with our experimental 
implementation. 

Example 8.8. Let TX = V{xl + y^,2xoXi + 'iylyi) = V{fi,f2) C C the tangent 

bundle of X = V{x'^ + y^) C C^. The annihilator of in D{S) and the Bernstein-Sato 

polynomial of TX using the previous approach can be computed with the SINGULAR 

commands Sannf sVar and bf ctVarAnn. 

LIB "bfunVar.lib"; 

ring R = 0, (xO,xl,yO,yl) ,Dp; 

ideal F = x0''2+y0~3, 2*x0*xl+3*y0~2*yl ; 

bf CtVarAnn (F) ; 

The output is lengthy, hence we supress it. We obtain an ideal called LD with 15 generators 
and the Bernstein- Sato polynomial for TX, which looks as follows 

hTx{s) = {s + l)\s + \)\s -t- \)\s + \){s + + I). 

Analogously, one can consider the tangent bundle o{V{x'^+y^). In this case the Bernstein 
polynomial has degree 42 and it equals 

4r(-ir('4r(-ir(-^)(-3)(^-A)(-^) 

2 
3 



12 
2 

15 



15/ 
5 

12. 



) 



5 

12 



The result was not able to be obtained using the elimination ordering given in section 
8.2. Instead the following monomial ordering (in SINGULAR format) has been taken: 
"(a(l,l),a(0,0,2,l,l,2),(dp(6),rp) ". 
Note that if TX is the variety defined by / = (/i, /2), then &/(s) has always negative 
roots. As this examples shows the same is not true for bz{s), due to the codimension. 



Example 8.9. Let Z be the algebraic variety defined by / 
x\x2)- Then 



{x\ — X2X3,X2 — XiX3,xl ■ 



This is actually Example 5.10 in Shibuta (2008) and it corresponds to the space of 

monomial ciirvc SpccC[T-^T'*,T'"^] C C^. Note that the 5-function coincides with the 
one that appears in Shibuta (2008), since we are computing hz{s) instead of 6/(s). 
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Example 8.10. Let Z be the so-called Hirzebruch-Jung singularity of type (5,2). It is 
a cyclic quotient singularity and can be seen as the algebraic variety associated with the 
ideal {z\ — Z2Z4, Z2Z3 — Z1Z4,, Z2 — ziZz) C C[2:i, z^, Z3, Z4]. Then 

bz{s) = {s+ins+his+l){s+h. 



Example 8.11. As an intractible example we would like to mention the following one. 
Let Z be the cyclic quotient singularity of type (6; 1, 2, 3). It can also be seen as the toric 
variety associated with the matrix 

/'642031 0^ 



A := 



12 3 10 
112 



yuuuuiizy 

The corresponding ideal can be taken to have 9 generators in C[zi, . . . , 2:7]. We are 
not yet able to compute the Bernstein-Sato polynomial even trying several elimination 
orderings. 



9. Conclusion and Future Work 

From our recent articles there follow some important conclusions. 

(1) In Levandovskyy and Martm-Morales (2008) we proved, that for the computation 
of Bernstein-Sato polynomial of a hypersurface, the homogenization, used in the 
method of Oaku and Takayama is superfluous. This, however, does not apply to 
the situation of computing a 6- function with respect to weights. By several steps in 
Levandovskyy and Martin- Morales (2008); Andres et al. (2009) and in this article 
we have shown, that the method by Briangon and Maisonobe can be seen as natural 
refinement of the method of Oaku and Takayama. 

(2) From our investigations it follows, that we cannot in general prove, that either 
initial-based method or annihilator-based one for the computation of Bernstein- 
Sato polynomial is definitely more efficient than the other one. Instead, on numer- 
ous examples we see that roughly the domain of better performance of initial-based 
method includes hyperplane arrangements, while for other singularities annihilator- 
based method scores distinctly better. It is important to continue these investiga- 
tions and derive more classes of polynomials, when possible, in order to use this 
information in attempts to estimate at least the practical complexity of f-module 
computations. 

(3) Also for the syzygy-driven method to compute Ann£)[s](/*) we do not have yet a 
proof of its superiority over the method by Briangon and Maisonobe. Due to the 
reasons we explain in this paper one could achieve such superiority. But on the 
other hand, there are examples, which show the contrary. Even if there are only a 
few examples of such kind, it is interesting to investigate this phenomenon deeper. 

(4) We pay so much attention to the algorithms for the case of a hypersurface due to 
many reasons. According to our proofs for the case of a variety, we use indeed the 
same technology. We expect to generalize all of enhancements we have described 
to the case of an afhne variety in a similar way we presented the gcineralization 
of algorithms for annihilator and Bernstein-Sato polynomial. It is very important 
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to stress, that working in the case of a variety is a priori much more; involved 
computationally, hence more attention on very effective algorithms need to be paid. 
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